PREPRINT UPR 0126MT 

June 1993 



UNIVERSITY of PENNSYLVANIA 
Department of Physics 
David Rittenhouse Laboratory 
Philadelphia PA 19104-6396 



From Skyrmions to NN phaseshifts 

Niels R. Walet 



Submitted to Physical Review C 



From Skyrmions to NN phaseshifts 

Niels R. Walet*'t 

Department of Physics, University of Pennsylvania, Philadelphia PA 19104-6396 

February 9, 2008 

Abstract 

We study the NN phaseshifts in a Hamiltonian obtained ft-om quantization of the 
collective modes in the Skyrme model. We show that a combination of an adiabatic and 
diabatic approximation gives a good NN force, with sufficient attraction to produce a 
bound deuteron. The description of the repulsive core appears to be the main cause for the 
remaining discrepancies between the Skyrme model force and phenomenology. Finally 
we discuss the possibility of finding non-strange dibaryon resonances in the J'^ = S"*" 
channel. 

1 Introduction 

In this paper we continue the study of the derivation of a nucleon-nucleon force from the 
Skyrme model. In previous work jl], ^ we have shown that the Skyrme model gives a good 
qualitative description of the intermediate and long range nuclear force. Here we study the 
phaseshifts obtained in the model. This requires new information about the repulsive core 
of the potential. Many more details of the problem, as well as some alternative approaches 
are reviewed in Refs. ^, ^ |5|. 

In our previous work we have discussed all the steps that are necessary to obtain a good 
result for the NN force. The starting assumption is that the Skyrme model a non-linear 
field theory of interacting pions, is related to QCD in the limit of a large number of colors 
(Nc). Even though this connection is uncertain, one can take the point of view that in the 
long wavelength limit only the lightest degree of freedom, the pion, plays a role. Thus one 
expects the model to be good for phenomena where only distances larger than the wavelength 
of the p meson, the next lightest particle, play a role. 

An interesting property of the Skyrme model is that there exists a topologically conserved 
current, usually associated with the baryon number B. The B = 1 solution, the Skyrmion or 
hedgehog, has been used to describe properties of hadrons to better than 30 % « 1/Ac, giving 
numerical support to the model. Once we take the Skyrme model to be a valid starting point 
for the discussion of single baryons, the next step is to study the interaction of two baryons. 
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To this end one studies a collective-coordinate manifold of B = 2 Skyrmion states, which is 
rich enough to describe the dynamics of NN states. Initial studies of the B = 2 system were 
based on the so-called product Ansatz, which is still used as a paradigm for understanding 
the structure of the manifold (see the discussion in Ref. 0). This Ansatz describes the 
superposition of two Skyrmions without any distortion. The advantage of this Ansatz is 
that it allows us to construct all states on the collective manifold in terms of a few invariant 
combination of six parameters. This in turn can be used to determine the collective potential 
and kinetic energy on all points of the manifold. To understand why the product Ansatz is 
not good enough for a realistic calculation of the A^A^ force, consider the potential energy in 
a little more detail. It is calculated as the difference between the energy of two interacting 
Skyrmions and twice the energy of a single Skyrmion. This number is typically of the order 
of tens of MeV's, but is obtained from the difference of two numbers of the order of 2 GeV. 
This is obviously sensitive to relatively minor flaws in the Ansatz. The product Ansatz has 
several flaws, that make it unlikely that it can lead to a correct description of the B = 2 
system. These are first of all that the product Ansatz is not symmetric under interchange 
of the two hedgehogs, secondly that it can not describe the cylindrically symmetric ground 
state of the B = 2 manifold, the donut ^, 10, 11]. 

In order to find improved solutions one can try to find better solutions on part of the 
collective manifold (the calculations are too demanding to solve for every state). To that 
end one picks out those sets of solutions that have a definite reflection symmetry (see, e.g., 
Ref. [|l2|). The Skyrmion configurations are now obtained numerically by imposition of a 
constraint on the separation of the two solutions, which leads to distorted states, finally cul- 
minating in a state of toroidal symmetry, the donut. This method is based on the assumption 
that these solutions are stationary against deformations that break the reflection symmetry. 
From preliminary results of a calculation performing a local RPA on the Atiyah-Manton 
Ansatz, this assumption appears to be largely correct [^]. As stated above the disadvantage 
of this approach is that one obtains only partial information about the collective Hamilto- 
nian. Fortunately the information one can extract coincides with the dominant part of the 
interaction as calculated in the product Ansatz. 

After performing all this analysis, we still end up with a classical collective Hamiltonian 
that describes the large limit, where the nucleon and its excited states are degenerate. 
Since we believe that the quantum world has = 2> we cannot requantize the collective 
Hamiltonian as it stands. Instead we requantize the Hamiltonian in a space of nucleons 
and A's. We use an algebraic technique to perform this reduction [14, |l^. This leads to 
a coupled channels problem, and not yet to a NN force. In a flnal step we invoked the 
Born-Oppenheimer approximation in our previous work to obtain the A^A^ adiabatic nuclear 
force, which was shown to be similar to a phenomenological force, such as the Reid soft-core 
potential |16|. 

In this paper we shall study the phaseshifts calculated in the full {N — A) case and for the 
adiabatic reduction. This will give us information about the adiabatic reduction, as well as 
about the quality of the intermediate coupled channels problem. We will have to go slightly 
beyond the adiabatic approximation in order to obtain information about the short-range 
repulsive core. This last part of the potential must be included in a calculation of NN 
phaseshifts. 
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Tes paper is organized as follows. In Sec. ^ we study our model for the collective Hamil- 
tonian and the reduction to a NN force. In the next section, Sec. D (^), we discuss the 
phaseshift calculation for the ^Sq and ^Si—^Di channels. Some technical details about the 
calculations can be found in the appendix. In Sec. § we discuss the possibility of finding 
resonances near A production thresholds. Finally, in Sec. we draw some conclusions. 

2 The model and the NN reduction 
2.1 The Skyrme model 

The Skyrme model |l^, ^] is a non-linear field theory that can be realized in terms of an 
SU(2)-valued matrix field U, with Lagrangian density 

C = ^TT[d^U{x)d>'U^{x)] + ^Tvp^d^U, U^d.Uf + ^[Tr(C/ + C/t) - 4]. (1) 

The model is covariant, as well as invariant under global SU(2)-rotations that are identified 
with the isospin symmetry. As was discovered by Skyrme the model has a topologically 
conserved quantum number, which is identified as the baryon number B. The U field is 
interpreted as a combination of a scalar a field and an isovector pion field, U = {a + iT-n)/ f.,^. 
The a field is not an independent physical field due to the unitarity constraint on U. 

The standard time-independent solution to the classical field equations for B = 1 is the 
defensive hedgehog, where the pion field points radially outward, 

Uiif) = exp{iT-ff{r)). (2) 

The baryon number of this state is given by = (/(O) — f{oo))/TT = 1. This solution breaks 
translational invariance, as well as the 0(4) spin- isospin symmetry. If we perform a global 
SU(2) isorotation on the state, 

Uiif^A) = A^Ui{r}A, (3) 

we obtain a state of the same energy. 

In the B = 2 system we will frequently use the product Ansatz, which is also used as a 
model to understand the more complicated numerical solutions. This Ansatz makes use of 
the fact that the product of two B = 1 solutions has baryon number two. The most general 
Ansatz we can construct from two hedgehogs consists of the product of two separated and 
rotated hedgehogs, 

U2{r\RAB) = A^Ui{r- R/2)AB^Ui{r + R/2)B 
= U2{r\RCD) 

= D^C^/^^Ui{f-R/2)CUi{f+R/2)C^/^^D. (4) 

In the last line of (Q) we have introduced the matrix D, that describes the rigid isorotation 
of the whole system, as well as a relative isorotation C. When R is very large changing C or 
D does not change the energy of the solution. For smaller R, D still generates a zero-mode 
(corresponding to broken isospin symmetry), but the energy will depend on C. Again, the 
energy is also invariant under spatial rotation, due to the conservation of angular momentum 
J = L + S. 
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2.2 Modelling the force 

As discussed in our previous work, the introduction of collective coordinates leads to an effec- 
tive Hamiltonian for the relative motion and relative orientation (we shall use the word rela- 



tive orientation as a synonym for the relative isorotation) . In the calculations of Refs. 1 1£ , 2C ] 
only a small part of the collective surface was calculated: Using reflection symmetries one 
studies three one-dimensional "lines" on the collective manifold where the radial coordinate 
R changes, but the orientation (analogous to the matrix C in the product Ansatz (^) re- 
mains fixed. This allows one to extract limited information about the collective Hamiltonian, 
through an expansion in invariants to first order. Fortunately the product Ansatz, where the 



complete expansion can be calculated [15|, allows one to show that higher order terms are 
small. 

The resulting classical Hamiltonian, obtained after making some crude approximations 
for the kinetic terms, still treats radial motion on the same footing as isorotational motion, 
which separates the nucleon from the A and higher I = J resonances. As is well known the 
appearance of unphysical resonances with higher / = J is due to the close relation of the 
Skyrme model with the string theory that describes the large Nc limit of QCD. In order to 
make progress towards the physical world, we need to re-impose the restriction Nc = 3. The 
way we proceed (see for an alternative but equivalent approach) is through a requantization 
of the spin-isospin degrees of freedom of the individual Skyrmions using an S'C/(4) interacting 
boson model for those degrees of freedom. The important step [|^, |l^ is to identify the boson 
number in these models with the number of colors. The representation space then consist of 
twenty states: the nucleon and A with all possible spin-isospin projections. 

The requantized Hamiltonian can be given in the following form 

Ti"^ 1 

H = Eo + — {R-^dnR^dR + L^/R^) + — {if + ll + Sf + Sl) + vi{R) + V2 {R)W + {R)Z{R) . 

(5) 

In the potential energy (the last three terms in Eq. (^)) the functions Vi are taken from the 
calculation of the interaction energy of two Skyrmions. The operators W and Z are similar 
to the spin-spin and tensor operators: 

W = T-T^jNl 

Z = T^Tl^^mRj-6i,]/Nl (6) 

Here a and f3 label two different sets of bosons, used to realize the m(4) algebras, and T is a 
one-body operator with spin and isospin 1. In order to simplify the algebra, we introduce a 
spherical tensor notation, which allows for the use of the standard Racah calculus. We then 
need the spin-isospin doubly reduced matrix elements of (the spherical tensor form of ) T. 
These are given by 

(iv|||r|||iv) = -10, (7) 

(A|||r|||A) = -20, (8) 
(A|||r|||iV) = -8^2. (9) 



4 



Table 1: Parameters in the fit to the potential Vi, as evaluated by The unit for r is fm, 
all energies are given in MeV. 



i 


1 


2 


3 


a\ 




1091.7 


1618.7 


a1 




-3156.4 


-4687.1 


b\ 


-5810.5 


8964.3 


13268.7 


hi 


22937.2 


117.47 


478.23 




-33370.8 






bf 


17532.7 







In the kinetic energy we have made the approximation that the mass is a constant, equal to 
twice the relevant reduced mass (so M still is an operator in isospace). We assume that both 
the mass M and the moment of inertia A are independent of R, and that A is independent of 
the orientation - certainly a gross simplification. Some preliminary work has been done on 
lifting those approximations |21|, but not enough is known at this instant to make a better 



approximation . 

As in our previous work we extracted the functions Vi from the work in Ref. Instead 
of using a table of numbers, we chose to make a fit to the results. One reason is that we 
shall need to push the calculation to short distances. We choose a form of powers of r times 
exponentials of one- and two-pion range: 

Vi{r) = {aj + af/r) exp{-r'm.„) + {h]r + b}) exp(-r2m^), z = 2, 3, (10) 
vi{r) = {b\r + bl + h\/r + b\/r'^)e^Y>{-r2m^). (11) 

The relevant parameters are listed in Table |l[ 

Note that the repulsive core is not very well determined from the standard calculations 
[^. This is due in part to the uncertainty in the definition of the coordinate R for small 
separations [R < 1.3 fm). 



2.3 Calculation of the A^A^ potential 

The calculation of an A^A^ potential from Eq. (^) can be split into two parts. One is the 
adiabatic calculation, valid when the A^A^ force is weak compared to the total energy of the 
system (from a numerical calculation one finds that, roughly, R > 1 fm). As discussed in 
Ref. ||2j, this is based on the separation of energy scales: For large Nq, one can distinguish 
two energy scales or reciprocally two time scales. The slower time scale is associated with 
the motion in the collective manifold, i.e., R and the orientation. The other time scale 
corresponds to the almost instantaneous response of the pion field to changes in R and the 
relative isospin orientation. For large A'^c we cannot separate the time scales for the two sets of 
adiabatic modes, as can be seen in the highly correlated doughnut. In the donut state orbital 
modes and isorotational modes are intrinsically linked, leading to the interesting structure 
of this state. 
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For Nq equal to three the situation changes. The R motion is typically much slower than 
the rotational motion which leads to the separation of the nucleon and A states. We thus 
have three energy scales, that of the pion field, of the — A separation and of the R motion. 
We can now calculate a Born-Oppenheimer potential for the i?-motion, which constitutes the 
slowest degree of freedom |]|, |] . 

This adiabatic - Born-Oppenheimer - approximation fails at shorter distances (i? < 1 fm) , 
where the R dependent potential becomes strongly repulsive. Numerical experimentation 
showed that for the calculation of the phaseshifts a diabatic prescription seems to work well. 
The adiabatic potential energy curves exhibit (relatively) narrow avoided crossing in the 
interior region. If we follow the adiabatic curve in the interior region, we obtain a potential 
of rather poor quality. If we follow the unmixed curve through the crossing, as if no crossing 
took place, we obtain much better results. This corresponds to a diabatic approximation, 
where we replace avoided crossings by real crossings. 



3 Calculation of phaseshifts 

We calculate the phaseshifts by integrating the Schrodinger equation from r = rmin where 
we impose hard-core boundary conditions (the forces are so repulsive at short distances, that 
this is no limitation), to r = rmax where we match to spherical Bessel functions describing 
the asymptotic behavior of the wavefunctions (this corresponds to the approximation that 
the potential vanishes beyond this point, but keeps all the slowly decaying centrifugal forces.) 
As discussed in the appendix this leads directly to an evaluation for the S matrix, which in 
its turn leads to the phaseshifts. Since we need to impose hard core boundary conditions at 
1^ = fmin, this point had better be inside the repulsive core of the potential. The choice of this 
interior radius is studied in fig. ||, where we show the calculated phaseshift in the ^Sq channel 
for center-of-mass energy of 1 MeV. The ordinate shows the radius rmax where we perform 
the matching, and the abscissa gives the calculated phaseshift. As can clearly be seen in the 
upper part of the figure, a hard-core radius of 0.7 fm gives the wrong phaseshifts. To be on 
the safe side, we have always used a radius of 0.3 fm in our calculations. Strengthened by the 
apparent robustness of our method, even close to threshold, we now study the phaseshifts 
in the ^5*0 channel as a function of the center-of-mass kinetic energy (the zero of energy is 
the A^A^ threshold). In fig. |2| we show the phaseshifts obtained from the full A^A coupled 
channels problem and compare it to the phaseshifts for the adiabatic A^A^ potential. It is 
important to note how close these results are. This implies that the adiabatic potential is a 
good approximation to the underlying more complicated problem. For good measure we have 



also added the Arndt et al phenomenological phaseshifts |2^, 23 1 obtained from the SAID 
database for np scattering in the same figure. Clearly we do not have enough attraction in 
this channel, even though the trend is surprisingly good. Actually we can obtain enough 



attraction by increasing all parameters in the potential function Vi in Eq. (11), except 64, by 
only 6 %! (Keeping 64 fixed means that we do not strengthen the most repulsive part of the 
interaction beyond its current value; if we increase this parameter we find it hard to obtain 
agreement.) As can be seen from Fig. ^, this gives enough attraction near threshold. The 
remaining difference between our model potential and the phenomenological phaseshifts is 
now mainly due to the description of the repulsive core. This is a part of the Skyrme model 
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Figure 1: The phaseshifts in the ^Sq channel {E = 1 MeV) as a function of the matching 
radius where we match to the Bessel function solutions of the free problem. We compare 
calculations for three different hard core radii. 
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Figure 2: Phaseshifts in the ^So channel as a function of center of mass energy. We compare 
the results with and without explicit A's to the phenomenological (VPI) phaseshifts. 
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Figure 3: Phaseshifts in the So channel after rescahng the potential as discussed in the text, 
we compare with the phenomenological phaseshifts and the unsealed results. 



that we have not yet analyzed. 

Of course we can perform the same analysis for the J = 1,T = 1 channel. We analyze 
the two-by-two S matrix in terms of the standard Stapp parameters 6± and e as [p^] 

( cos2ee^'^- i sin 2e e*(^++^-) \ , . 

~ zsin2e e^('5++^-) cos2ee2^'5+ )' ^ > 

Again the results from the full potential and the adiabatic approximation are very close. So 
close that we have not plotted them in fig. ^. As can be seen the results compare favorably 
with the phenomenological analysis, even though we have less attraction. To our surprise, 
however, the phase-shift goes to 180° at threshold. This shows that our potential supports 
a bound state. Further analysis shows that this occurs very close to threshold, the state 
is bound by only 0.053 MeV, reflecting the fact that the attraction in the current channel 
is too weak. Again, a simple rescaling by 6 % of all parameters but 64 gives a binding 
energy much closer to experiment, Eb = 1.77 MeV. If we improved on the calculation of the 
repulsive core (which appears to be too strong in the current calculation) one might even 
expect better results, since both the phaseshifts and the binding energy of the deuteron are 
obtained through the delicate balance between attraction and repulsion. This makes it even 
more surprising that we obtain results so close to experiment. 



4 Behavior near A production thresholds 



the T matrix 



In many of the phaseshift analyses by the VPI group (e.g., Refs. [22, |' 
exhibits a resonance very near the A production threshold. As better data became available 
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Figure 4: The mixing parameter ei and the phaseshifts in the deuteron {^Si—^Di) channel 
as a function of center-of-mass energy. We compare the phase shifts of the coupled A^A 
problem (labeled full) with the phenomenological phase shifts of the VPI group. 
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Figure 5: The A — A potential in J = 3^, L = channel, compared to the NN potential in 
the J = 1^, L = channel. Neither potential includes effects from the Born-Oppenheimer 
approximation . 



some of these poles have disappeared, but currently a few channels still seem to exhibit 
resonances [^]. The abovementioned analysis has not yet been pushed to near the double 
A production region. This region is of some interest, however. 

Theoretically it has been argued that strong attraction in the J = 3^ channel is "in- 



evitable" [g5|. Recently Wang et al studied this question using an approach to baryon 
interactions based on nonrelativistic quark dynamics. Using their "quark delocalization" 
scheme and color screening, they find a A — A potential in the S-wave, J'^ = 3"*", T = chan- 
nel that is attractive and 400 MeV deep. Our model can be used to study the same question, 
since we have explicit A's in our calculations. Using the same non-relativistic Hamiltonian 
as used in the previous section we can discuss the problem of the A — A potential. 

In Fig. |5| we show our result for the A — A potential in the J'^ = 3"*", T = channel as 
a function of separation. We see that there is indeed an attractive potential in this channel, 
but that it is no where near the unusually large potential reported by Wang et al. In this 
case our potential is obtained without channel mixing (since we are studying excited states 
mixing will actually decrease the binding). If we perform a coupled channels bound-state 
calculation in the A — A subspace only we find a bound state with a binding energy of about 
5 MeV. 

We have solved for the nucleon-nucleon phaseshifts below the A — A threshold in the 
r = 0, J'^ = 3"*" channel, where the only two open channels are the A^A^ D and G waves. Our 
results are only qualitative, since we use a non-relativistic coupled-channel analysis of our 
potential, including the closed channels. (Note that Wang et al also used a non-relativistic 
model.) As can be seen in Fig. ^ we find a sharp resonance in the D-wave, 5 MeV below the 
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Figure 6: Phaseshifts in the NN J = 3^ channels just below the A — A threshold (at 600 
MeV in our calculation). 

A — A threshold. In order to see how robust such a resonance is we include the width of the 
A in our calculations. For a resonance close to threshold one can make the approximation 
that the threshold energy is complex, with an imaginary part corresponding to the width of 
the A. Using this approximation we find that the resonance is extremely sensitive to such a 
change. Even a width of 1 MeV is enough to destroy the resonant behavior. 

How relevant are these results? We find very reasonable (if somewhat weak) attraction 
in our NN calculation. It is thus likely that the relatively small attraction in the A — A 
channel is not too far from phenomenology. This would rule out an observable non-strange 
dibaryon once we take into account the width of the A. 

5 Conclusions 

In this paper we have shown that the Skyrme model can be used to obtain a NN interaction 
that is reasonable in terms of the phase shifts it produces. The calculations reported in this 
paper give phaseshifts that are surprisingly close to phenomcnological ones. With minor 
modifications we can obtain suflacient attraction to get both the deuteron binding energy 
and the attraction in the ^Sq channel close to the experimental value. The weakest point of 
the potential is actually both the description and the size of the repulsive core. This should 
not surprise us since there is a great ambiguity in the definition of the radial coordinate 
at small distances, which leads to a great ambiguity in the size of the repulsive potentials. 
More importantly the Skyrme model itself should not be trusted at such distances. This last 
problem can probably be solved by adding more vector mesons to the model (the minimal 
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extension is probably the inclusion of the /?, u; and Ai). The first problem may prove harder 
to solve, even though we have recently made some progress on this problem using the Atiyah- 
Manton Ansatz in combination with techniques of large-amplitude collective motion. This 
same approach may be used to shed some light on the change of inertial parameters with 
distance and orientation, which leads to a modification of the kinetic energy part of the 
Hamiltonian |13]. In order to understand the behavior of the inertial parameters it would be 



useful to understand their behavior in the product Ansatz first. This may lead to clues about 
which terms are important, and which are not. A comprehensive study of these parameters 



in the product Ansatz is currently in progress |l21 |. 
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A Calculation of phaseshifts 

To calculate the phaseshifts we express the Hamiltonian in a coupled channels matrix form, 
(after re-expressing the equations in terms of Ui = rTpi) 

Here we require that Vij{oo) = 0, and thus the threshold for channel i is given by Ef^'^'^^ . 
The integer Lj is the orbital angular momentum of the given channel In order to solve the 
scattering problem, we impose hard-core boundary conditions at a small radius rmin. We 
also assume that the potentials Vij are negligibly small for a radius rmaxj so that we can 
approximate solutions by free solutions for r larger than this value. 

We integrate outwards from rmin; given the derivative at this point. For convenience we 
introduce Vi{r) = drUi{r). The differential equation then reads 



drUi{r) = Vi{r), 
drVi{r) = —2mi ^ 



E - £;thres Li{Li + 1) 



2m,- 



Ui{r)-Vi,u,ir)] . (14) 



We solve this by an explicit method (i.e., one that only uses function values at previous 
values of r to determine what happens at the next value). If we keep the initial derivative 
o = ^''(rjnin) as an explicit parameter, this can be written in matrix form 

(15) 

where the matrix M depends on the method used (we use a fourth order Runge Kutta). Due 
to the hard core boundary conditions we thus only need two of the sub-matrices, 

u = M^'^'a, V = M^'^'a. (16) 
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At this point we match to the spherical Bessel functions that describe the asymptotic behavior 
of the problem, i.e., the eigenvalues of a Hamiltonian consisting only of kinetic terms and 
centrifugal forces. 



Ui{r) = aikirh^l]{kir) + I3ikirh}i'.{kir). 



(17) 



Here the channel momentum ki = y2mi{E — Ef^'^'-^^), with /cj in the positive half-plane if 
E < _£;*'hresi_ -yy^g j^^^ solve foF Q!j and Pi by matching at r = rmaxj (It is convenient to use 
the notation h^{kir) = hrh^l^ (kir)) 



(Xi = 



(2), 



= M+a, 



A = 



-/i?'(fciW)fciM«" - h>l>{kir^)kiM 



= M_a. 



(18) 



We have to impose the condition that there are no outgoing waves (which blow up exponen- 
tially) in the closed channels, i.e., for those values of i where ki is complex. Let us label the 
closed channels by c and the remaining open channels by o. The condition that A is in the 
closed channels can easily be implemented: 







can be shown to imply, by inversion of the matrix, that 



a = 



1 T\/rcofn/roo h/Toc 1 n/rco\ — li:io 



We can thus relate a° to 



Mr - 



Ml 



-M" 



^3 



Here we have used the fact that our choice of phases leads directly to the S matrix. 



(19) 



(20) 



(21) 
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